EVOLUTION OPERATORS FOR QUANTUM CHAINS 
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Abstract. Discrete-time evolution operators in integrable quantum lattice models are sometimes 
more fundamental objects then Hamiltonians. In this paper we study an evolution operator for the 
one-dimensional integrable g-deformed Bose gas with XXZ-type impurities and find its spectrum. 
Evolution operators give a new interpretation of known integrable systems, for instance our system 
describes apparently a simplest laser with a clear resonance peak in the spectrum. 

Existence of a complete set of commuting operators is the basic principle of quantum integrability. 
However, often the commuting set does not provide a distinguished operator deserving the title of 
the Hamiltonian for a physical system. In particular, this is the common feature of quantum models 
obtained by quantization of classical equations of motion in wholly discrete space-time. Equations 
of motion define a discrete time translations for classical variables, a map A(r,a) — > A(t + l,cr). 
Corresponding translation for quantum observables is produced by an evolution operator, 

(1) A(Y + l,cr) = UA^trJUt , UU f = 1. 

If the time unity interval coincides with the unity spacing, then there is no a small parameter expansion 
of U defining a lattice Hamiltonian. The evolution operator becomes the main object of the lattice 
quantum mechanics. 

Discrete-time evolution operators for quantum chains were considered in many papers. For instance, 
the study of the spectrum of evolution operator was in the focus of the quantum Liouville theory pp . In 
this paper we discuss another example of quantum evolution system: the one-dimensional q-deformed 
Bose gas j2j with XXZ-type impurities. The simple test "what evolution operator is doing in the space 
of states" will give a new view on well known models. 

Quantum mechanics begins with the algebra of observables. In our case the local algebra of observ- 
ables is the g-oscillator j^] with generators B, B^, N: 

(2) BBt = 1 - g 2N + 2 , BtB = l-<? 2N , 0<«<1. 
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Below we use the Fock space representation of J2J), the Fock vacuum is defined by B|0) = 0. The 
^-oscillators are the algebra of observables for q-deformed Bose gas. In addition, it is very convenient 
to use the g-oscillator representation for XXZ model |S] . In this paper we consider the combined 
chain of three types of ^-oscillators, 

(3) {B jini Bt 3> , N 3> } : j = 1,2,3, n = 0, 1, . . . , N - 1 , 

where index 3 „ stands for the g-Bose gas, indices x,r» aud i.n stand for oscillator representation of XXZ 
model in n th site of the combined chain. The quantum Lax operators are 

~N 3 , n uB f 



-q UU'z.n 



(4) M n (u) 
for the g-Bose gas (condition q < 1 is equivalent to the attractive potential), and 

(5) L n (u) = 



uq 



for the XXZ model. L-operator (JSJ) has the center 

(6) s„ = i(Ni, n + N 2 , n ) 

with the half-integer spectrum, matrix (JSJ on the subspace s„ = s n is identical to the usual spin-s„ 
L-operator for U q {sl2). Note, on the subspace s n = operator Q is just unity matrix, L n (u)\ Srz =o = 
(u + e' e ) 1. Extra parameter e in JSJ) is equivalent to an exponential shift of the spectral parameter u. 
The transfer matrix of the model is defined by 

(7) t{u) = TracerMo(M)io(w)Mi(M)Li(u)M 2 (M)i2(u)---Mjv_i(u)iiv-iK 

Complete set of commutative operators is given by the decomposition of t(u) and by the set of charges 
s„ ©. Subjects of our interest are sub-sectors of the extended Hilbert space corresponding to specified 
eigenvalues of s„. In particular, if all s„ = 0, this is just the pure one-dimensional g-Bose gas with 
attractive potential. If So = 1/2 and all other s„ ~ 0, this is the Bose gas with the single impurity. If 
some of s„ = 1/2 and all the other s„ = (our final case), this is the q-Bose gas with some density of 
impurities. 

Our main object, the discrete-time evolution operator, is defined by 



(8) 



UM n (u)L n (u) = L n {u)M n+l {u)V 



Evidently, transfer matrix 10 commutes with the evolution operator and therefore gives the integrals 
of motion. One can rewrite (jHJ as the map: 

U ^."B^IT 1 = q N ^B\, n - e ie g Nl -Bt 2: „B 3 ,„ , 

(9) U Bt^U" 1 = B+ liTl B+3,n + qe' e q N ^ + ^B^, n , 

U g N2 -+ 1 Bt 3i „U- 1 = g Nl ."+ 1 Bt 3l n+i - e ie g N3 -+ 1 B 1 , n+1 Bt 2)n+1 . 
Corresponding relations for UB^ n U _1 are the conjugation of 10: U is the unitary operator if q and 
£ are real. Equations J^J provide in addition Us„ = s„U and UK = KU where 

(10) K = (N 2 ,„ + N 3 ,„) . 

n 

Note, the map @ is exactly the tetrahedral map from All two-dimensional models have hidden 
three-dimensional structure. 

Take up now the interpretation of the evolution operator. In what follows, we use 



(11) q = e 
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We start with the lowest occupation numbers over the total Fock vacuum |0). Let U|0) = |0), and 
we consider next the subspace of the Hilbert space with Sq = K = 1. Let us see, what the evolution 
operator makes with the state 

(12) = Bt 2i0 |Q> . 
Applying © once, we get 

(13) U - e i£ "^ |^) + |Vo) ® \b ) , 
where we use the notation 

(14) |Vo) ® |6*> = Bt li0 Bt 3jfc |0> . 
What is happening after r steps of discrete evolution: 

(15) U r |^) = e^^ T + e-< ie -^ |^o) ® ( ^ e^-^^lbk) J . 

\fc=o / 

Suppose, r is big enough, but 

(16) 7T <C 1 , er ~ some integer multiple of 2tt . 

Then equation Ijl5|) describes the decay the quasi-stable state with the energy e and width 7. The 
second term in (|15fl is a state with the bosonic wave function ^(r, k) — e( l6_7 ^ r_fc ) - the speed-of-light 
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right-moving wave radiated by excited state. Due to (|16J) . the oscillating term in the wave function 
dominates, and therefore the wave package in the right hand side of 1)15(1 is the state with the average 
energy and momentum e. Note, we may talk about the energy of the radiated state only if r in (|15fl 
is big enough, this exactly corresponds to the Heisenberg uncertainty principle. The one-step distance 
of the lattice may be considered as the ultra-microscopic scale of the space-time, while big r in (|15|l 
corresponds to a microscopic scale. 

Thus, the discrete-time evolution operator provides an alternative terminology. We will call site n 
of the spin chain with s n = 1/2 as the atom with two energy states: the ground state \ij)o) with zero 
energy and the excited state with the energy e and width 7. The regime for 7 and e follows from 
®: 

(17) 7 < e < 2tt . 

The spin chain sites with s n = do not provide any effect, these sites are just empty. The Bose field 
Bt 3 n is identified with the electro-magnetic field, the radiated state in i|15|) is the photon. The whole 
system is evidently a kind of chiral laser (photons are radiated to the right only). Note, the impurities 
destroy the translation invariance of the system. Scattered and radiated photons may have momenta, 
but the atoms are rigidly fixed and the solid media takes up the recoil momenta of photons. 

Consider next eigenstates and eigenvalues of the evolution operator. We consider the periodical 
boundary conditions - a "toroidal chiral laser" . The eigenstates have the structure of the superposition 
of stationary bosonic waves scattering on the atoms (the coordinate Bethe Ansatz). For instance, the 
one-boson and one-atom eigenstate such that U|* p ) = e lp |\I>j,), is 

(18) \%) = qeip [ _J J Bt 2 , o |0) + Bt 1)0 Bt p |0) , 

where 

(19) Bt.^e-^Bt^ 

is the creation operator of the photon. The quantization equation for the momenta of one-photon 
states follows from the periodical boundary conditions, 

pip _ np' s 

(20) e' N P 6 . g6 = 1 . 

Since the evolution operator describes only the right-moving photons with dispersion relation E = p, 
the momentum p in (|19I2U|) must be positive. 



In general, the spectrum of U is given by 



k ( K 



(21) 



U = Y[ui = exp i \ y^Pi 



where Ui = e' Pi are the roots of the Bethe Ansatz equations 



(22) 



q 2s "Ui — e' e 



n 



q 1 u i - quj 
qui - q~ x Uj 



i,j = 1, ...,K . 



n 



Here N is the length of the chain, periodical boundary conditions are taken into account, and K is 
the eigenvalue of (|10fl (the number of photons in the laser). In general, s n are arbitrary eigenvalues of 
s„, but here we consider s n = 0, 1/2 only. Equations l|22|l are literally the Bethe Ansatz equations for 
transfer matrix (0; relation (|21|l is the main analytical result of this paper. 

The eigenstates corresponding to (|22|l may be interpreted in terms of the photon creation operators 
in momentum space (|19|) . If all pi are real and different, the photonic counterpart of the eigenstate is 
a slightly modified state B^B^ ...'B ji PK |0). In addition to real Pi, the Bethe Ansatz equations (|22|l 
have solutions with complex p t (strings of the Bethe Ansatz equations) , for instance p\ ~ p + \"f and 
p 2 ~ p - tj. The bosonic counterpart of such eigenstate is e ~'P( k + k )-i\k-k |gt 3 ^.B^.^'IO). This is 



the bound state for big 7, but in our regime l|17f) it may be understood as properly 7-regularized state 
Btp|0). Analogously, the highest bound states correspond to the highest occupation numbers B^lO). 
The important point is that we have to solve the Bethe Ansatz equations implying the right-moving 
condition Re(p;) > 0. 

Turn finally to the thermodynamical limit N,K — > 00. Let the number of two-states atoms be 
N' < N, their density n = N'/N is not zero at the limit. Structure of the infrared part of spectrum 
is the subject of separate investigation. The lowest energy levels correspond to Bethe Ansatz strings, 
the 7-regularization solves apparently the infrared catastrophe problem. 

In this paper we consider the "optical wave band" p ~ e. Assume that the Bethe roots Ui = e' Pi , 
lm(pi) = 0, form a dense distribution Q{pi)~ l = N(pi+± — p%) near pi ~ e. Assume in addition that 
we can neglect in this region densities of photons with higher occupation numbers (complex pi). The 
standard Bethe Ansatz technique provides the Hulthen integral equation for the density Q(p): 



(23) 




where 
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Solving (|23|l in the optical wave band, we neglect the contribution of the radio band since — S(p,p') — > 

dp 

when \p — p'\ ^7. The solution is 

27 

what is evidently the laser resonance peak over the white noise. Accounting of radio band and complex 
Pi may modify the white noise term. The resonance term has exactly the AAX-spin chain structure, so 
that the number of resonant photons in the stationary state is half of the number of atoms. Additionally 
pumped energy does not increase the density of resonant photons, it dissipates into the white noise. 
Acknowledgements. I would like to thank M. Bortz, J. de Gier and V. Mangazeev for fruitful 
discussions. 
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